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Abstract 

Some consequences of spatio-temporal symmetry for the deterministic 
decomposition of complex light fields into factorized components are con- 
sidered. This enables to reveal interrelations between spatial and temporal 
coherence properties of wave. An estimation of average number of the de- 
composition terms is obtained in the case of statistical ensemble of light 
pulses. 
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1. INTRODUCTION 

Modal description of light coherence, being a multidimensional generalisation 
of the well-known Karhunen-Loeve expansion, was first introduced in optics by 
Gamo|l|. In short, the essence of the approach lies in the fact that any correlation 
function of a field E{r, t) — in particular the transverse beam coherence rs(r, r') 
— can be expressed as a superposition of factorized components 

Ts{v, v') = [ dtEir, t)E*{r', t) = ^ M„^„(r) S^ir'), (1) 

where each term of the sum represents a completely coherent partial wave Sn{r). 
The decomposition basic functions and modal energies m„ are eigenvectors and 
eigenvalues of Fredholm's integral equation 

u^Sn{r)= j £r'Vs{r,v')£^(j') (2) 
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and, since the kernel of (2) is Hermitian, the set of functions Sn{r) is orthonormal 

J d'^r£nir)£^{r) = 6n,m- (3) 

In (1) the transverse correlation function r5(r, r') is determined as a time 
average (over the pulse duration or time of registration) and, in this sense, it 
is a deterministic characteristic of any particular wave. In the case of statistical 
ensemble of similar pulses the time averaging can be replaced by statistical one or 
both types of averaging may be combined. Formally it has no effect on relations 
(1) - (3), but, as it will be seen later, changes their physical meaning. 

The modal expansion (1) - (3), as well as its modification for space-frequency 
domain]^, ^, is broadly used in coherence theory as a convenient tool for estima- 
tion of entropy and informational capacity of light beams[0], as the best method 
for modelling of coherence properties of a complex light wave with a finite set of 
simple mutually incoherent waves[|], and so on. These relations constitute the 
mathematical basis for proof of various types of uncertainty inequalities]^]. 

The next principal step in the development of the modal formalism was made 
by Pasmanik and Sidorovich[^]. They demonstrated the spatio-temporal symme- 
try for decomposition (1) - (3), that leads to some important relations between 
spatial and temporal coherence characteristics of light waves. The discussion of 
them both in deterministic form and under application of the ensemble averaging 
is a main goal of the present paper. 



2. DUAL APPROACH OVERVIEW 

Let us suppose that for some wavefield E{r,t), where r is a two-dimensional 
radius-vector at a plane z = const, the modal expansion (1) - (3) is known. So 
far the solutions of the integral equation (2) form a complete functional basis 
(when including the functions, corresponding to zero eigenvalues^), one can 
define a set of projections of the initial field onto this basis 

v/^e„(t) = J d'rE{r,t)£:{r), (4) 

that, in turn, allows to build up a modal representation of the field itself 

E(r,t) = ^V^^„(r)e„(t). (5) 

n 

The most important point here is that for pulses of finite total energy 

U = J d'^rj dt\E{r,t)f < oo 

the projections e„(t) also constitute the complete orthonormal set of basic func- 
tions 

J dt en{t)e*^{t) = 6n,m- (6) 
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Last relation can be proved by direct substitution of definition (4) into (6) and 
accounting (3). It is just this mutual orthogonality of temporal functions e„(t), 
that leads to absence of any interference between different terms of spatial basis. 

Another approach to evaluation of temporal basis (4) lies in use of dual integral 
equation 

UnCnit) = J dt'TT{t,t')en{t'), (7) 

where rT{t,t') is a global temporal correlation function of the field E{r,t) 

TTit,t') = [ d\ Eiv,t) E*{r,t') = J2ur, er.it) elit'). (8) 

In contrast to the standard definitional (with averaging over the time or en- 
semble of pulses) the averaging procedure in (8) is carried out over the beam 
cross-section. Hence, the function rr(t,t') expresses the overall correlation the 
wavefront patterns and is closely related to the degree of similarity^, ^ H{t,t') 
of the wavefield for consequent time moments 
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H{t,t') 



J d^rE{r,t)E*{r,t') j d^r |^(r,t)|') [j d^r \E{r,t') 



One can easily see that equations (2) and (7) make up two equivalent dual 
variants for evaluation the decomposition (5). Both equations have identical 
spectra of eigenvalues and for complete description of modal structure of field 
one needs to know only one set of basic functions Sn{r) or e„(t). The second can 
be immediately determined through projection (4) or via its dual equivalent 

T^Snir) = J dtE{r,t)e*St). (9) 

The last variant (6) - (9) has an advantage of dealing with 1-D task. One 
more exact consequence of the dual formalism is that two functions rs(r, r') 
and Vxit^t') in the case of no degeneracy (all Un are different) allow one to 
completely reconstruct the form of field. Under degeneration (e. g. when Ui = Uj) 
the ambiguity arises from the fact that two different wave structures ei{t)£i{r) + 
ej{t)£j{r) and ei{t)£j{v) + ej{t)£i{r) produce the same correlation functions (1), 
(8). 

Strictly speaking, the modal structure of the field does not remain constant 
under the wave propagation, but mode mixing is comparatively low for quasi- 
monochromatic beams with small divergence [^]. That is why the deterministic 
dual decomposition (5) is inherently aimed to description of laser pulses and has 
been first applied in nonlinear optics0, where the partial coherence just means 
a high complexity of interacting waves. 

In practice the complete modal description can be fulfilled only for very few 



classes of models||T^, |TT|, what is, first of all, related with intricacy of integral 
equations (2), (7) solving. Therefore those consequences of the method are taking 
the special significance, for which one does not need to know the exact basic 
functions e„(t), £^n(r). 
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3. EFFECTIVE NUMBER OF MODES 



So far as the mode number n in general cannot be univalently associated with 
any other parameter of partial wave (except its energy), the only natural way 
to restore distribution of without solving (2), (7), is evaluation of nonlinear 
A;-order moments of modal spectra J^n'^n- can be done with use of iterated 
kernels theorem via sequential integration of functions r5'(r, r') or rT{t,t'). The 
moment of 1st order has a trivial meaning of total field energy 

U = J2un= f d'rTsir,r) = f dtTrit^t). 

As it shown in Ref. 12, the higher moments determine a probability distri- 
bution of wave amplitude under conditions of wave mixing at a strong scatterer. 
The most important characteristic of field structure is an effective number of 
terms in decompositions (1), (5), (8)[|l], ||, [l^, which is expressed through 2nd 
order moment 

Neff = {UnUnf I (En M^) 

= U^ I {jd^Jd'r' \rs{r,r')\')=U' / {j dt J dt' Mt,t')l')- ^^^^ 

The value of Ngff specifies the ability of the total field to produce interference 
effects between two arbitrary separate points of the beam cross-section p, and 
changes from unity for spatial coherent one-mode wave to infinity for completely 
incoherent field. 

Two equivalent forms of (10) reflect real interconnection between spatial and 
temporal parameters of a beam. If one determines effective area of beam cross- 
section — Sf.jf, area of spatial coherence — cTc, pulse duration — T^jj and corre- 
lation time — Tc in the form 

Scff = (/ d'rTsM^' I (/ rfVr|(r,r)) , (ll.a) 

T,g = (I dtTTit,t)y I {j dtY%{t,t)\^ , (11.6) 
ac=[j d\ I d'p \Ts{r + p/2,r-p/2)f^ j [j dVr|(r,r)), (ll.c) 
'JdtJdT\Ts{t + T/2,t-T/2)f'^ I (^j dtTl{t,t)y (ll.rf) 



then relation (10) takes the form of equality for spatial and temporal degrees of 
freedom of wavefield 

^ = ^. (12) 
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It means that number of coherence zones per beam cross-section is equal to 
number of different spatial patterns over the pulse duration. 

Three of introduced in (11) parameters — T^ff, Sejj, ac — have quite tradi- 
tional meaning [0 and need no special remarks. The averaged over beam cross- 
section coherence time Tc describes time of global changing of field structure or, 
in other words, characteristic width hy t — t' of the degree of similarity t') of 
spatial wave patterns. Definitions (10), (11) have no sensitivity to overall phase 
modulation of the field 

E{Y,t) ^ E(r,t)exp(i0(r) +#(t)) 

and, therefore, value of Tc can rather significantly differ from local correlation 
time, which is defined in signal theory. 



4. APPLICATION OF STATISTICAL AVER- 
AGING 

Till this point the basic formalism has dealt with a wave- field -E(r, t) as with the 
deterministic one. At the same time classical coherence theory usually operates 
with radiation characteristics, averaged over ensemble of similar fields, because in 
the majority of cases the individual pulse parameters are not of interest. Hence, 
the natural question arises — how can such stochastic hypotheses influence on 
the results of previous analysis? 

As a first step let us consider what statistical averaging gives at the stage of the 
basic integral equations (2), (7) formulation. Just as in the standard approach fll], 
one can substitute the kernels r5(r, r') and Trit, t') with their averages (r5(r, r')) 
and {rT{t,t')). However it is evident, that transversal correlation function (1), 
averaged over time interval only, contains much more information about spatial 
structure of a beam, than similar value of (r5(r, r')) does. The approximate 
equality can take place only in the limit of infinite pulse duration and under 
quasiergodicity of the ensemble. Exactly the same with appropriate changing 
of words can be stated about temporal correlation function rT{t,t'). It is easy 
to see that such a lack of information breaks the main property of the present 
formalism — its spatio-temporal duality. 

As a result of this procedure the interpretation of the relations (1) - (4) and 
(6) - (9) must be changed. All other conclusions of the above sections (with 
exception of equality (12), which disappears) remain valid if taking into account 
that now we talk about two different and complementary means for describing of 
coherence of random pulse ensemble (but not for a particular wave). For transver- 
sal coherence this will be nonstationary variant of Gamo's treatment [jl| and for 
temporal correlation function it will have the form of modified Karhunen-Loeve 
expansion ||14|| with double averaging — over ensemble and beam cross-section. 
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In the subsequent discussion we shall utilise the fact that for each decomposition 
one can introduce its own number of degrees of freedom (10) — spatial Ns and 
temporal Nt 

^ (^') / (/ / I(r5(r,r'))r) = ^^lim ^iVeif, (13.a) 

Nt = {U') I (/ dt I dt' |(r^(t,t'))l') = j,}fi^N,ff- (13-&) 

Now it is clear that in order to preserve spatio-temporal duality of the formal- 
ism, the ensemble averaging should be applied at the later stages of consideration. 
As an example of such approach let us estimate the number of terms in the modal 
decomposition of a mean light pulse from the ensemble. The simplest and the 
most popular type of field statistics is Gaussian. In this case the coarse estimation 
can be done by averaging of integrals in definition of N(,jj (10), that, accounting 
the splitting of higher correlations and (13), gives a very simple formula 

On the basis of general reasons one can formulate some more requirements, 
which a priori should be satisfied by any admissible dependence Ngff(Ns, Nt). 
Thus, the function Neff{Ns, Nt) must be symmetrical about permutation of its 
arguments because of the dual status of spatial and temporal degrees of freedom 



N,ff{Ns, Nt) = N,ff{NT, Ng). (15) 



The value of Ngff must be a non- decreasing function of its arguments, that with 
(13) leads to conclusion 

'Kff < Ns, Nt 

and in asymptotics Ns — const • Nt ^ 1 the average number of modes will be 
linear with respect to any of arguments, in particular. 



NeffiNs = Nt = N :^ 1) (X N. 

Ensemble with only one degree of freedom in any of the basic subspaces corre- 
sponds to coherent (in terms of (5)) field 



NeffiNs = 1, Nt) = N,ff{Ns, TVt = 1) = 1- 

It is easy to see that estimation (14) obeys all above requirements but the last 
one, i. e. it poorly describes the region of small numbers of degrees of freedom (it 
is the consequence of approximate way of averaging (Neff)). The situation can 
be improved by taking into account the fiuctuations not only of the denominator, 
but also of the numerator (energy of light pulses) of expression (10) 

U^) = {U)'{l + e{Ns,NT)), 
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J (fr j £r' j dt j dt'\{E{T,t)E*{v\t'))\ 
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Correction e must satisfy the condition (15) and have the order of magnitude 
e oc l/(A^s'A^r), that can be confirmed by consideration cross-spectrally pure 
light [|I^, when correlation function factorizes. Hence, the refined estimation of 
the number of modes in the mean pulse may be written as following 

NsNt + 1 



Neff = ■ (16) 



For the first time estimation like (16) was given without a proof in paper[16 
for a system of several identical, statistically independent emitters with drifting 
phase. At limit Ns, Nt ^ 1 appropriate formula from Ref. 16 converts to 
(14). One can point out some more cases, which asymptotically lead to the same 
dependence. All this allows to say that area of applicability of relation (14) as 
estimation of N^ff is much wider than above assumptions. 

In order to illustrate the consequences from relation (14) we can consider a 
very vivid example of the ensemble of Schell-model fields ||T0|, 



{E{r, t)E*{r', t')) = ^/(r, t)I{r', t')j{r -r',t- t'). (17) 

One of the possible interpretations of model (17) corresponds to illustrative situ- 
ation when a fast shutter cuts off a pulse of radiation from a primary steady-state 
uniform partially coherent source. Then just beyond the shutter the degree of 
coherence 7(r,t) is specified by statistical parameters of the source only (say, 
with cr, r being an area and a time of correlation, respectively), while I{r,t) is 
(within a factor) a deterministic function of the shutter transmittance {S, T — 
the shutter aperture area and the time it is opened). On substituting (17) into 
(13) and accounting (11), (14) one can assure that the effective area of coherence 
(in the mean pulse) depends not only on spatial parameters, but as well on ratio 
between temporal characteristics T/r of the primary source and the shutter. And 
vice versa, the lifetime of a particular wavefront structure in the mean pulse is 
also determined by the ratio of S/a. It explains the significance of the discussed 
modal formalism for the nonlinear optical and laser beam problems. 



5. DISCUSSION 

In conclusion it is worth to point out the resemblance of the considered modal 
technique with bi-orthogonal decompositions used in other branches of physics 
— e. g. turbulence theory^ and pattern recognitionp^. Such tie is based on the 
common concept of complex process representation. By this analogy, the spatial 
partial coherence may be described as a sequence of more or less similar frames 
(instant wavefield structures) replacing each other. From this viewpoint coherent 
modes specify the feature basis of wavefronts evolution. 
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According to the general concept the apphcation of global ensemble averaging 
procedure is efficient (it gives results with comparatively small relative variance) 
when the number of modes in the mean pulse is high. Nevertheless, there are 
situations where under small N^ff the number of the ensemble degrees of freedom 
in one of the subspaces is much more than it [Nt ^ Nf.ff or Ns 3> N^jj). In this 
case the statistical averaging over the corresponding complex substructure of the 
field may be useful. 

How it is seen from (5), each single mode in deterministic decomposition 
produces factorized correlation function, i. e. corresponds to cross-spectral pure 
light But if we go to the whole multimode field, the spectral purity vanishes. 
Moreover, for the statistical ensemble even one-mode field will, in general, not be 
cross-spectral pure. 

Besides the discussed manifestations of spatio-temporal symmetry for modal 
decomposition, there exists a wide class of uncertainty relations]^, where it must 
also appear. They should have the form of inequalities bounding modal charac- 
teristics with such parameters of wave as angular divergence and spectral band- 
width. This statement leans against the fact that proof of uncertainty relations 
for correlation functions does not depend on the type of averaging used. 



References 

[1] H. Gamo, "Matrix treatment of partial coherence". Progress in Optics, Ed. 
E. Wolf, 3, pp. 187-332, North-Holland, Amsterdam, 1964 

[2] E. Wolf, "New theory of partial coherence in the space-frequency domain. 
Part I: Spectra and cross-spectra of steady-state sources", J. Opt. Soc. Am., 
72, pp. 343-351, 1982. 

[3] E. Wolf, "New theory of partial coherence in the space-frequency domain. 
Part II: Steady-state fields and higher-order correlations", J. Opt. Soc. Am., 
A3, pp. 76-85, 1986. 

[4] P. De Santis, F. Gori, G. Guattary, C. Palma, "Synthesis of partially coher- 
ent fields", J. Opt. Soc. Am., A 3, pp. 1258-1262, 1986. 

[5] M. J. Bastiaans, "Uncertainty principle for partially coherent light" , J. Opt. 
Soc. Am., 73, pp. 251-255, 1983. 

[6] G. A. Pasmanik, V. G. Sidorovich, "Interrelation between coherent proper- 
ties of light beams and space-time structure", Radiophys. & Quantum Elec- 
tron., 23, pp. 809-814, 1980. 

[7] G. Berkooz, P. Holmes, J. L. Lumley, "The proper orthogonal decomposition 
in the analysis of turbulent flows", Ann. Rev. Fluid Mech., 25, pp. 539-575, 
1993. 



8 



[8] A. M. Lazaruk, "Limits of applicability of the Raman-Nath approximation 
in problems of radiation self-diffraction", Optics & Spectroscopy (USA), 53, 
pp. 633-636, 1982. 

[9] B. Ya. Zel'dovich, N. F. Pilipetsky, V. G. Shkunov, Principles of phase 
conjugation, Springer- Verlag, Berlin, 1985. 

[10] A. Starikov, E. Wolf, "Coherent-mode representation of Gaussian Schell- 
model sources and of their radiation fields", J. Opt. Soc. Am., 72, pp. 923- 
928, 1982. 

[11] F. Gori, G. Guattary, C. Padovani, "Modal expansion for Jo-correlated 
Schell-model sources". Opt. Comm., 64, pp. 311-316, 1987. 

[12] A. A. Leshchev, "Method for measuring the number of spatially coherent 
modes in optical radiation". Optics & Spectroscopy (USA), 55, pp. 599-600, 
1983 

[13] A. Starikov, "Effective number of degrees of freedom of partially coherent 
sources", J. Opt. Soc. Am., 72, pp. 1538-1544, 1982. 

[14] Fu K. S. Sequential methods in pattern recognition and machine learning. 
Academic Press, New York - London, 1968. 

[15] J. W. Goodman, Statistical Optics, A Wiley-Interscience Publication, New 
York, 1985. 

[16] M. V. Vasil'ev, V. Yu. Venedictov, A. A. Leshchev, P. M. Semenov, 
V. G. Sidorovich, O. V. Solodyankin, "Reduction of speckle-structure con- 
trast in image under laser illumination". Optics k, Spectroscopy (USA), 70, 
pp. 6-7, 1991 

[17] J. Deschampus, D. Courjon, J. Bulabois, "Gaussian Schell-modcl sources: an 
example and some perspectives", J. Opt. Soc. Am., 73, pp. 256-261, 1983 



9 



